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In this paper we propose a dark matter model and study aspects of its phenomenol- 
ogy. Our model is based on a new dark matter sector with a U{\)' gauge symmetry 
plus a discrete symmetry added to the Standard Model of particle physics. The 
new fields of the dark matter sector have no hadronic charges and couple only to 
leptons. Our model can not only give rise to the observed neutrino mass hierarchy, 
but can also generate the baryon number asymmetry via non-thermal leptogenesis. 
The breaking of the new U(l)' symmetry produces cosmic strings. The dark matter 
particles are produced non-thermally from cosmic string loop decay which allows one 
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to obtain sufficiently large annihilation cross sections to explain the observed cosmic 
ray positron and electron fluxes recently measured by the PAMELA, ATIC, PPB- 
BETS, Fermi-LAT, and HESS experiments while maintaining the required overall 
dark matter energy density. The high velocity of the dark matter particles from 
cosmic string loop decay leads to a low phase space density and thus to a dark mat- 
ter profile with a constant density core in contrast to what happens in a scenario 
with thermally produced cold dark matter where the density keeps rising towards 
the center. As a result, the flux of 7 rays radiated from the final leptonic states 
of dark matter annihilation from the Galactic center is suppressed and satisfies the 
constraints from the HESS 7-ray observations. 

PACS numbers: 14.60.Pq, 95.35. +d 



I. INTRODUCTION 



There is strong evidence for the existence of a substantial amount of cold dark matter 
(CDM). The leading CDM candidates are weakly interacting massive particles (WIMPs), 
for example, the lightest neutralino in supersymmetric models with R parity. With a small 
cosmo logical constant, the CDM scenario is consistent with both the observations of the 
large scale structure of the Universe (scales much larger than IMpc) and the fluctuations of 
the cosmic microwave background 

However, the collisionless CDM scenario predicts too much power on small scales, such 



as a large excess o 



dwarf galaxies 



: dwarf galaxies ation of dark matter (DM) in 

6] and in large galaxies Q]. To solve this problem, two of us with their 
collaborators proposed a scenario based on non-thermal production of WIMPs, which can 
be relativistic when generated. The WIMPs' comoving free-streaming scales could be as 
large as or possibly even larger than 0.1 Mpc. Then, the density fluctuations on scales less 
than the free-streaming scale would be suppressed [8|. Thus, the discrepancies between the 
observations of DM halos on sub-galactic scales and the predictions of the standard WIMP 
DM picture could be resolved. 

Recently, the ATIC [9] and PPB-BETS [lo| collaborations have reported measurements 
of the cosmic ray electron/positron spectrum at energies of up to ~ 1 TeV. The data shows 
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an obvious excess over the expected background for energies in the ranges ~ 300 — 800 GeV 
and ~ 500 — 800 GeV, respectively. At the same time, the PAMELA collaboration also 
released their first cosmic-ray measurements of the positron fraction H] and the p/p ratio 
fl^ |. The positron fraction (but not the antiproton to proton ratio) shows a significant 
excess for energies above 10 GeV up to ~ 100 GeV, compared to the background predicted 
by conventional cosmic-ray propagation models. This result is consistent with previous 



measurements by HEAT 



13] and AMS 



Q. 



Very recently, the Fermi-LAT collaboration has released data on the measurement of the 
electron spectrum from 20 GeV to 1 TeV [l^], and the HESS collaboration has published 
electron spectrum data from 340 GeV to 700 GeV 16(, complementing their earlier mea- 
surements at 700 GeV to 5 TeV 17] . The Fermi-LAT measured spectrum agrees with ATIC 
below 300 GeV; however, it does not exhibit the special features at large energy. There have 
already been some discussions on the implications for DM physics obtained by combining 
the Fermi-LAT, HESS and PAMELA results Q. 

The ATIC, PPB-BETS and PAMELA results indicate the existence of a new source of 
primary electrons and positrons, while the hadronic processes are suppressed. It is well 
known that DM annihilation can be a possible origin for primary cosmic rays [3] which 
could account for the ATIC, PPB-BETS and PAMELA data simultaneously, as discussed 



first in 1 201 and also in 1211 (see 



22j for a list of references) 1 . However, the fact that the 



p/p ratio does not show an excess gives strong constraints on DM models if they are to 
explain the data. In particular, it is very difficult to use well-known DM candidates like the 
neutralino to explain the ATIC and PAMELA data simultaneously 24J since they would 
also yield an excess of antiprotons. Therefore, if the observed electron/positron or positron 
excesses indeed arise from DM annihilation, it seems to us that there may exist special 
connections between the DM sector and lepton physics 25j (see also 26l. 1271] ). 

In this paper, we propose a DM model and study its implications for DM detection. We 
fit our model to two different combinations of the experiment data: one set of data from 
the ATIC, PPB-BETS and PAMELA experiments; the other from the Fermi-LAT, HESS, 
and PAMELA experiments. Our results show that our model can naturally explain the e ± 



1 Note, however, that there are also astrophysical (see e.g. [23[) or other particle physics (see e.g. [22j]) 
explanations. 
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excesses while at the same time solving the small scale problems of the standard ACDM 
model via non-thermal DM production. For a single Majorana DM particle, its annihilation 
cross section has s wave suppression. Thus, we consider two degenerate Majorana DM 
particles. We add a new DM sector with a {7(1)' gauge symmetry and introduce an additional 
discrete symmetry to the Standard Model (SM). The DM particles are stable due to the 
discrete symmetry. During the {7(1)' gauge symmetry breaking phase transition a network 
of cosmic strings is generated. The decay of cosmic string loops is a channel for producing 
a non-thermal distribution of DM. This non-thermal distribution allows for DM masses 
and annihilation cross sections large enough to explain the cosmic ray anomalies while 
simultaneously remaining consistent with the observed DM energy density. In addition, the 
observed neutrino masses and mixings can be explained via the seesaw mechanism, and the 
baryon number asymmetry can be generated via non-thermal leptogenesis 28] . 

It has been recently recognized that a large annihilation cross section of DM particles 
into leptons to account for the cosmic ray anomalies will induce a large flux of 7 rays from 
the Galactic Center (GC) [29] or from the centers of dwarf galaxies 30]. The predicted 7 
ray fluxes based on the NFW profile for the standard CDM scenario have been shown to 
be in slight conflict with the current observations of HESS [si|. However, in our model the 
DM particles are produced non-thermally, so the high velocity of the DM particles will lower 



32]. 



the phase space density of DM and lead to a DM profile with a constant density core 
Therefore our model with non-thermally produced DM on one hand gives rise to a large 
annihilation cross section to account for the positron/electron excess observed locally while 
on the other hand it suppresses the DM density at the GC and leads to a low flux of 7 ray 
radiation. 

Our paper is organized as follows: in Section II, we describe in detail the model and 
the production mechanism of the DM particles. In Section III we study aspects of the 
phenomenology of the model, including studies of some constraints on the model parameters 
from particle physics experiments, implications for the PAMELA, ATIC, PPB-BETS, Fermi- 
LAT, and HESS results, and also the 7-ray radiation from the GC. Section IV contains the 
discussion and conclusions. 
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II. THE DARK MATTER MODEL 
A. The Dark Matter Sector 

The DM model we propose consists of adding a new "DM sector" to the Standard Model. 
The new particles have only leptonic charges and are uncharged under color. This ensures 
that the DM particles annihilate preferentially into leptons. To ensure the existence of a 
stable DM particle, the new sector is endowed with a discrete symmetry which plays a role 
similar to that of R-parity in supersymmetric models. The lightest particles which are odd 
under the Z 2 symmetry which we introduce are the candidate DM particles. 

In our convention, we denote the right-handed leptons and Higgs doublet as e R (l, — 1) 
and H(2, — \) = (H°,H~) T , respectively, where their SU(2) L x U(l)y quantum numbers 
are given in parenthesis. 

We consider the generalized Standard Model with an additional U(l)' gauge symmetry 
broken at an intermediate scale. In particular, all the SM fermions and Higgs fields are 
uncharged under this U(l)' gauge symmetry. To break the U(l)' gauge symmetry, we 
introduce a SM singlet Higgs field S with U{\)' charge —2. Moreover, we introduce four 
SM singlet chiral fermions xi, Xii N\, and N 2 , a SM singlet scalar field E and a SM doublet 
scalar field H' with SU(2) L x U(1) Y quantum numbers (1, —1) and (2, |), respectively. The 
U(l)' charges for Xi and H' are 1, while the U(l)' charges for iVj and E are — 1. Thus, our 
model is anomaly free. To have stable DM candidates, we introduce a Z 2 symmetry. Under 
this Z 2 symmetry, only the particles Xi an d E are odd while all the other particles are even. 
The x particles will be the DM candidates, whereas the chiral fermions iVj will play the role 
of right-handed neutrinos. 

The relevant part of the most general renormalizable Lagrangian consistent with the new 
symmetries is 

-£= \m^S + + \m\,H*H> + ±(S*S)* + ^Ef + ^ H'f 

+^-(S^S)(E^E) + ^(&E)(H*H') + ^-(S^S)(H'^H') + ^-(S^S)(H^H) 

Zi _ _ 

+ypx!Xj + yZ&NiNj + y^L i H'N j + H.c.) . (1) 
As we will discuss in the following subsection, the vacuum expectation value (VEV) for 



6 



S is around 10 9 GeV. Then, the couplings A3, A 5 and A 6 should be very small - about 1CT 12 
- in order for the model to be consistent with the expected value of the SM Higgs. This fine- 
tuning problem could be solved naturally if we were to consider a supersymmetric model. 
Moreover, in order to explain the recent cosmic ray data, the Yukawa couplings should be 
around 10 -6 . This would generate a DM mass around 1 TeV. Such small Yukawa couplings 



ylg can be explained via the Froggat-Nielsen mechanism |33j which will not be studied here. 

To explain the neutrino masses and mixings via the "seesaw mechanism" , we require that 
the VEV of H' be about 0.1 GeV if y 1 ^ ~ 1 and y l J ~ 1. In this case, the lightest active 
neutrino is massless since we only have two right-handed neutrinos iVj. In addition, in our 
i7(l)' model, the Higgs field forming the strings is also the Higgs field which gives masses 
to the right-handed neutrinos. There are right-handed neutrinos trapped as transverse 
zero modes in the core of the strings. When cosmic string loops decay, they release these 
neutrinos. This is an out-of-equilibrium process. The released neutrinos acquire heavy 
Majorana masses and decay into massless leptons and electroweak Higgs particles to produce 
a lepton asymmetry, which is converted into a baryon number asymmetry via sphaleron 



transitions 28J. Thus, we can explain the baryon number asymmetry via non-thermal 
leptogenesis. 

In this paper, we consider two degenerate Majorana DM candidates Xi an d X2 since 
the annihilation cross section for a single Majorana DM particle is too small to explain 



the recent cosmic ray experiments 25| . For simplicity, we assume that the Lagrangian is 



invariant under xi ^ X2- Thus, we have 

yl = y ':, y« - y* . (2) 



To make sure that we have two degenerate Majorana DM candidates xi an d X2, we choose 
y^ = 0, and assume m x < rag. 



B. Non-Thermal Dark Matter Production via Cosmic Strings 

We assume that the {7(1)' gauge symmetry is broken by the VEV of the scalar field S. 
To be specific, we take the potential of S to be 

V(S) = ~\(\S\ 2 - V 2 ) 2 , (3) 
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where A is the self-interaction coupling constant. The VEV of S hence is (S) = rj with m| = 
At? 2 . Due to finite temperature effects, the symmetry is unbroken at high temperatures. 
During the cooling of the very early universe, a symmetry breaking phase transition takes 
place at a temperature T c with 

T c ~ V\r). (4) 

During this phase transition, inevitably a network of local cosmic strings will be formed. 
These strings are topologically non-trivial field configurations formed by the Higgs field S 
and the U(l)' gauge field A. The mass per unit length of the strings is given by \i = rf. 

During the phase transition, a network of strings forms, consisting of both infinite strings 
and cosmic string loops. After the transition, the infinite string network coarsens and more 
loops form from the intercommuting of infinite strings. Cosmic string loops loose their 
energy by emitting gravitational radiation. When the radius of a loop becomes of the order 
of the string width w ~ A _1 / 2 r/ _1 , the loop releases its final energy into a burst of A and 
S particles 2 . Those particles subsequently decay into DM particles, with branching ratios 
e and e' . For simplicity we assume that all the final string energy goes into A particles. A 
single decaying cosmic string loop thus releases 

N ~ 2ii\- 1 e (5) 

DM particles which we take to have a monochromatic distribution with energy E ~ ^P, the 
energy of an S-quantum in the broken phase. In our model, we assume that the masses for 
A, S and Ni are roughly the same, so we have e = 1. 

Given the symmetry we have imposed, the number densities of xi an d X2 are equal. Thus, 
the number density Udm °f DM particles, the sum of the number densities of Xi an d X2, is 

n D M = n xi + n X2 = 2n Xl = 2n X2 . (6) 

If the S and A quanta were in thermal equilibrium before the phase transition, then the 
string network is formed with a microscopic correlation length £(t c ) (where t c is the time at 
which the phase transition takes place). The correlation length gives the mean curvature 



2 We are not considering here DM production from cosmic string cusp annihilation since the efficiency of 
this mechanism may be much smaller than the upper estimate established in [jyj, as discussed e.g. in 
351 ] . DM production from cusp annihilation has been considered in 36 1 . 
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radius and mean separation of the strings. As discussed in [37| (see also the reviews [38| ) . 
the initial correlation length is 

fa) ~ A" V 1 . (7) 
After string formation, there is a time interval during which the dynamics of the strings 
is friction-dominated. In this period, the correlation length increases faster than the Hubble 
radius because loop intercommutation is very efficient. As was discussed e.g. in [39[, the 
correlation length scale £(i) in the friction epoch scales as 

3 

t \ 2 



m = m i -j . (8) 

The friction epoch continues until £(t) becomes comparable to the Hubble radius t. After 
this point, the string network follows a "scaling solution" with £(t) ~ t. This scaling solution 
continues to the present time. 

The loss of energy from the network of long strings with correlation length £(£) is pre- 
dominantly due to the production of cosmic string loops. The number density of cosmic 



string loops created per unit of time is given by [38|, |39( : 



dn A df 

IE = ■ < 9 > 

where v is a constant of order 1. We are interested in loops decaying below the temperature 
T x when the DM particles fall out of thermal equilibrium (loops decaying earlier will produce 
DM particles which simply thermalize). We denote the corresponding time by t x . 

The DM number density released from t x till today is obtained by 8] summing up the 
contributions of all decaying loops. Each loop yields a number N of DM particles. We track 
the loops decaying at some time t in terms of the time tf when that loop was created. Since 
the loop density decreases sharply as a function of time, it is the loops which decay right 
after t x which dominate the integral. For the values of Gfi which we are interested in, it 
turns out that loops decaying around t x were created in the friction epoch, and the loop 
number density is determined by inserting ([8]) into (Q. Changing the integration variable 
from t to £(£), we integrate the redshifted number density to obtain: 



t 



n n D T(to) = Nu J ^-J T 4 ^ , (10) 

where the subscript refers to parameters which are evaluated today. In the above, £p = 
£,(tp) where tp is the time at which cosmic string loops which are decaying at the time t x 
formed. 
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Now the 
at a rate 



oop's time-averaged radius (radius averaged over a period of oscillation) shrinks 

-JT = —FloopsGfi , (11) 

where Ti oops is a numerical factor ~ 10 — 20. Since loops form at time tp with an average 
radius 

R(t F ) ~ X^g^GuMffi, (12) 

where g* counts the number of massless degrees of freedom in the corresponding phase, they 
have shrunk to a point at the time 

t * AV^M. (13) 

Thus 

t F ~ x- 1/3 r l/2 rf oops M n hl (14) 

Now the entropy density is 

S = ^T\ (15) 

The time t and temperature T are related by 

-i Mpi 

t = 0.3^ 2 (T)-^, (16) 

where Mpi is the Planck mass. Thus using Eqs. (jHJ) and fflUj) . we find that the DM number 
density today released by decaying cosmic string loops is given by 

^nonth f\ 7K T' 4 

ynonth _ »DM _ ^L,.\ 3/2x^-2 3/2 „-5/2 w2 lX 

r DAf — s — ^ CZ/A 1 loops9*T c 9*T x 9*T F IV1 PlrpQl \ LI ) 

where the subscript on g* refers to the time when g* is evaluated. 
The DM relic abundance is related to Y x by: 

-U2 



Q x h m m x Y x s(t )p c (t ) h 



2.82 x 10 8 Y* ot (m x /GeV) , (18) 
where h is the Hubble parameter in units of 100kms _1 Mpc~ 1 , m x is the DM mass, and 

~ytot ~ytherm _i_ ynonth 

X ~ X ' X 

To give some concrete numbers, we choose the parameter values e = 1, u = 1, X = 0.5, 
r = 10, Mpi = 1.22 x 10 19 GeV and Q x h 2 = 0.11. In our model, we have #* T = 136, 
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TABLE I: The required T c values in units of GeV for various choices of a and j3 in the cases 
m x = 620 GeV, m x = 780 GeV, and m x = 1500 GeV, respectively. 



a 


1 


1 


2 


2 


5 


5 


P 


1 


0.5 


1 


0.5 


1 


0.5 


T c (m x = 620 GeV) 


7.7 x 10 9 


8.6 x 10 9 


4.8 x 10 9 


5.4 x 10 9 


2.6 x 10 9 


2.9 x 10 9 


T c (m x = 780 GeV) 


9.3 x 10 9 


1.0 x 10 10 


5.9 x 10 9 


6.6 x 10 9 


3.2 x 10 9 


3.6 x 10 9 


T c (m x = 1500 GeV) 


1.6 x 10 10 


1.8 x 10 10 


1.0 x 10 10 


1.1 x 10 10 


5.5 x 10 9 


6.2 x 10 9 


a 


10 


10 


15 


15 


20 


20 


P 


1 


0.5 


1 


0.5 


1 


0.5 


T c (m x = 620 GeV) 


1.7 x 10 9 


1.9 x 10 9 


1.3 x 10 9 


1.4 x 10 9 


1.0 x 10 9 


1.2 x 10 9 


T c (m x = 780 GeV) 


2.0 x 10 9 


2.2 x 10 9 


1.5 x 10 9 


1.7 x 10 9 


1.3 x 10 9 


1.4 x 10 9 


T c (m x = 1500 GeV) 


3.5 x 10 9 


3.9 x 10 9 


2.6 x 10 9 


3.0 x 10 9 


2.2 x 10 9 


2.4 x 10 9 



g* Tp = 128, and g* T ^ =128. We define the dimensionless ratios 

ynonth 

= (19) 

Demanding that we obtain a specific value of f3 for the above choices of the parameter values 
will fix T c via f|T8|) . For various values of a and ft, we present the resulting T c values for 
the cases m x = 620 GeV, m x = 780 GeV, and m x = 1500 GeV, respectively, in Table fl] In 
short, T c must be around 10 9 GeV if we want to generate enough DM density non-thermally 
via cosmic strings. 



III. PHENOMENOLOGY OF THE MODEL 

A. Constraints on the Model Parameters 

The coupling constants y % e between right-handed leptons and the DM sector are con- 
strained by experiments, and especially by the precise value of muon anomalous magnetic 
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moment g — 2. Assuming that the masses of x an d E are nearly degenerate, we obtain that 
the contribution to the muon anomalous magnetic moment from the new coupling is about 





1 m 2 



Sa ' ~ ' < 20) 

X 



4l|, 



The 2cr upper bound from the E821 Collaboration on 5a ^ is smaller than ~ 40 x 10 10 
from which we get for m x ~ 1 TeV, 

2/„<10. (21) 
For the electron anomalous magnetic momentum we assume the contribution from the dark 



sector is within the experimental error 



42| 



5a e < 7 x 10"" . (22) 

Then we get a upper limit on y e which is about 30. Therefore the constraints on the couplings 
of the model due to the heavy masses of the new particles are quite loose. 

Now we study the constraints from the experimental limits on lepton flavor violation 
(LFV) processes such as /i — > ej, r — > //(e)7 and so on. The branching ratios for the 



radiative LFV processes are given by 40J 



Br{e t - en) ~ a em m*/2 x (^g^) /I\ , (23) 

where Tj is the width of ej. Given the experimental constraint on the process /i — > we 
get 

Sr(/i _ e7 ) „ lo-s x (y e y^ < lO" 11 , (24) 
which gives that y e y^ < 0.03. For the process r — > n(e)j we have 

5r(r -> /i(e) 7 ) ~ 10~ 9 x (y, (e) y T ) 2 < lO' 7 , (25) 

which leads to the conclusion that y T y^( e ) J$ 10. Connecting the DM sector to the PAMELA 
and Fermi-LAT (or ATIC) results usually requires a large branching ratio into electron and 
positron pairs. From the LFV constraints shown above we conclude that it is possible to 
have a large branching ratio for the annihilation of the DM particles directly into e + e~ , or 
via A i+ /i". 



12 



B. Explanation for the Cosmic e Excesses 



In our model the DM sector only couples to the SM lepton sector. Therefore DM anni- 
hilates into leptons dominantly. Furthermore, since DM is produced non-thermally in our 
model the DM annihilation rates can be quite large with a sizable Yukawa coupling y\. Thus 
our model can naturally explain the cosmic excesses observed. 

Because the annihilation cross sections for X1X1 an d X2X2 to leptons are s wave suppressed, 
the dominant cross sections of X1X2 annihilating into charged leptons are given by [2J 



\vi\ 2 \vi\ 2 - 



32tt 



x In 



Sa Is (s — 4m 2 ) 



s (s — 4m 2 ) + 



s + 2m ~ — 2m% 



s + 2m| — 2m 2 — y s (s — 4m 2 ) 



s + 2m 2 ~ — 2m 2 + J s (s — 4m 2 ) 



x 



2\ 2 



2 (m| - ml) 



,(26) 



s + 2m~ — 2m 2 — y s (s — 4m 2 ) s + 2m~ — 2m 2 + y s (s — 4m 2 ) 

where v is the relative velocity between the two annihilating particles in their center of mass 
system. The overall factor 4 will be cancelled when we calculate the lepton fluxes, so, we 
will leave it in our discussions. Up to 0(v 2 ), the above cross section can be simplified as 25] 



a {j v ~ 



1287T 



lv£lV.I 2 



+ 



(2 + r) 2 L(2 + r) 2 (2 + r) 3 _ 



m: 



(27) 



where 



2 2 
m ~ — mz 

r= E „ x >0. 



(28) 



With f ~ 10 3 and r ~ 0, we obtain 



25] 



(a ijV ) < 4 x 1.2 x lO^Wsec" 1 (Z9°M') j^Vj 2 



m. 



(29) 



We emphasize that the Yukawa couplings should be smaller than y 47r for the perturbative 
analysis to be valid. 

In our model with non-thermal production of DM particles, we consider two separate fits 
to the ATIC /PPB-BETS /PAMELA and Fermi-LAT/HESS/PAMELA datasets. Firstly we 
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consider a numerical fit to the ATIC, PPB-BETS and PAMELA data 25j. In this case we 
assume the DM mass to be 620 GeV and that DM annihilates into electron/positron pairs 
predominantly, i.e., y\ ~ for i — 2, 3. In the second case we fit the Fermi-LAT, HESS 
and PAMELA data by taking the DM mass 1500 GeV and assuming that DM annihilates 
into pairs dominantly. Note that all lepton fluxes resulting from DM annihilation are 

proportional to n 2 cr ann for models with a single DM candidate \- Because n Xl = n X2 = n x /2 
in our model, the lepton fluxes are proportional to 

1 

4' 



(30) 



This will cancel the overall factor 4 in the above annihilation cross sections in Eqs. fl26l) and 

dm 





bkg 


bkg + 0.62 TeV(eV) 


bkg + 1 .5 TeV(u V 




ATIC i— •— i 


- PPB-BETS i — A — i 




HESS w—> 


t T l ~TT l HESS 


jw energy 


i — 1 — I — t-rij Fermi LAT i — * — i 














t\l i 


jl 



bkg - 

bkg + 0.62 TeV(eV) ■ 
bkg + 1 .5 TeV(iiV) ■ 
A MS v 
HEAT94+95 t 
HEATO0 
PAMELA > 




FIG. 1: Left: The e + + e spectrum including the contribution from DM annihilation compared 
with the observational data from ATIC [9], PPB-BETS [10], HESS Q, Q and Fermi-LAT Q. 
Right: The e + /(e~ + e + ) ratio including the contribution from DM annihilation as a function of 



3, 



43|] and PAMELA 



111 ]. Two sets of 



energy com pared with the data from AMS Q, HEAT 
fitting parameters are considered: in one model (Model I) the DM mass is 620 GeV with e + e~ 
being the main annihilation channel to fit the ATIC data, while in the other model (Model II) 
the DM mass is 1500 GeV and we assume that is the main annihilation channel to fit the 

Fermi-LAT data. 



In Fig. [T]we show that both cases can give a good fit to the data after considering the 



propagation of electrons and positrons in interstellar space 25| with the annihilation cross 
section 0.75 x lCE 23 cm 3 s -1 and 3.6 x 10 -23 cm 3 s _1 , respectively. The model parameters of 
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the two fits are given in Table HT1 For the first fit, we do not need the boost factor at all by 
choosing y\ = 2.6, which is still smaller than the upper limit V^7r for a valid perturbative 
theory. Moreover, choosing y\ = 3 in the second fit, we just need a small boost factor about 
10 which may be due to the clumps of the DM distribution [44|. Therefore, the results on 
the observed cosmic excesses can be explained naturally in our model. 



C. 7- Ray Radiation from the Galactic Center 

Since the explanations of the anomalous cosmic ray require a very large annihilation 
cross section to account for the observational results, this condition leads to a strong 7-ray 



radiation from the final lepton states. In particular, observations of the GC 29[ or the center 



of dwarf galaxies 30J have already led to constraints on the flux of the 7-ray radiation. 

The HESS observation of 7-rays from the GC [31J sets constraints on the Galactic DM 
profile. The NFW profile in the standard CDM scenario leads to too large a flux of 7-rays, 
thus conflicting with the HESS observation. On the other hand, if DM is produced non- 
thermally as suggested in Section II the DM profile will have a constant density core 32| so 
that the 7-ray radiation from the GC will be greatly suppressed. 

In our numerical studies, we consider the following two cases to constrain the DM profile: 

• Case I: we simply require that the 7-ray flux due to final state radiation (FSR) do not 
exceed the HESS observation. 

• Case II: we make a global fit to the HESS data by assuming an astrophysical source 
with power law spectrum plus an additional component from FSR resulting from DM 
annihilation. 

Let us consider a DM profile taking the form 

p(r) = - Ps , (31) 



where p s is the scale density and r s = r vir /c V i r (l — 7) is the scale radius, with r vir the virial 
radius of the halo 3 and c v i r the concentration parameter. In this work the concentration 



3 The virial radius is usually denned as the range inside which the average density of DM is some factor of 

n A /(i+2) 3 +n A 



the critical density p c , e.g., 187r 2 + 82a; — 39a; 2 with x — CIm(z) — 1 = —75 — TTTz^Tn ^ or a ACDM universe 



45] 
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FIG. 2: Upper: the FSR 7-ray fluxes from a region with \l\ < 0.8° and \b\ < 0.3° close to 
GC compared with the observational data from HESS [31]. The left panel compares the two 
models given in Table [II] directly with the data, while the right panel shows the combined fitting 
results using a power law astrophysical background together with the FSR contribution from DM 
annihilation at 95% (2<r) confidence level. Lower: constraints on the DM profile parameters 7 and 
c v i r due to the HESS observation of 7-ray radiation from the GC by assuming different final leptonic 
states. The left panel corresponds to the constraint Case I, while the right panel corresponds to 
Case II. The two curves in the right panel represent the la and 2a upper bounds respectively. 



parameter c v i r and shape parameter 7 are left free, and we normalize the local DM density to 
be 0.3 GeV cm -3 . Then the virial radius and total halo mass are solved to get self-consistent 
values. Given the density profile, the 7-ray flux along a specific direction can be written as 

<j)(E,^) = C x W{E) x J(V>) 

_„!«o x Wdiv 1 / A 



4tt 2ml dE p%R & J LOS 
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where the integral is taken along the line-of-sight, W(E) and J{ip) represent the particle 
physics factor and the astrophysical factor respectively. Thus, if the particle physics factor 
is fixed using the locally observed e + e~ fluxes, we can get constraints on the astrophysical 
factor, and hence the DM density profile, according to the 7-ray flux. For the emission from 
a diffuse region with solid angle AQ, we define the average astrophysical factor as 

Jmi = -^J J{ii>)dQ. (33) 

The constraints on the average astrophysical factor Jaq for the two models are gathered 
in Table [III in which JXn* shows the maximum J factor corresponding to Case I, while J/^ a 
corresponds to Case II, at the 68% (la) and 95% (2er) confidence levels. The 7-ray fluxes 
of the two cases are shown in the upper panels of Fig. [2j 

In the lower panels of Fig. [2] we show the iso- J^n lines in the 7 — c V { r plane for Case 
I (left) and Case II (right) respectively. In this figure we also show the mass condition of 
(1 - 2) x 10 12 M of the Milky Way halo. From Fig. El we can see that the NFW profile with 
7=1 (chosen based on N-body simulation in the standard CDM scenario) is constrained 
by the HESS data, if the observed cosmic e ± excesses are interpreted as DM annihilation. 
However, if DM is produced non-thermally the high velocity of the DM particle will make 
the DM behave like warm DM and lead to a flat DM profile which suppresses the 7-ray flux 
from the GC. 

TABLE II: Parameters of the two scenarios adopted to fit the ATIC /PPB-BETS /PAMELA or 
Fermi-LAT/HESS /PAMELA data. 





channel m x (GeV) (av}(W 23 cm 3 s 1 ) 


rmax 
J AQ 


Tier 
J AQ 


t2(T 


Model I 


e + e- 620 0.75 


300 


42 


97 


Model II 


H+jT 1500 3.6 


200 


81 


111 



IV. DISCUSSION AND CONCLUSIONS 

In this paper we have proposed a DM model and studied aspects of its phenomenol- 
ogy. We have shown that our model can simultaneously explain the cosmic ray anomalies 
recently measured by the ATIC, PPB-BETS and PAMELA experiments or by the Fermi- 
LAT, HESS and PAMELA experiments, resolve the small-scale structure problems of the 



17 



standard ACDM paradigm, explain the observed neutrino mass hierarchies, explain the 
baryon number asymmetry via non-thermal leptogenesis and suppress the 7 ray radiation 
from the GC. 

In this model, DM couples only to leptons. In direct detection experiments it would 
show as an "electromagnetic" event rather than a nuclear recoil. Experiments that re- 
ject electromagnetic events would thus be ignoring the signal. However, in the Fermi- 
LAT/HESS/PAMELA fits, the DM particle couples mainly to muons, and there being no 
muons in the target of direct detection experiments, no significant signal would be expected. 
In the ATIC /PPB-BETS /PAMELA fit, the DM couples predominantly to electrons; the 
electron recoil energy is of order m e v^, M ~ 0.1 eV, and it would be too small to be de- 
tectable in current devices. Alternatively, this energy could cause fluorescence [46], albeit 
the fluorescence cross section would be prohibitively small. Regarding the annual modula- 
tion signal observed by DAMA 47J, although this experiment accepts all recoil signals, an 
estimate of the electron scattering cross section shows that the present model predicts a cross 
section which is about 8 orders of magnitude smaller than ~ lpb required to account for the 
modulation [26]. Therefore we do not expect a signal in direct detection experiments if the 
DM model presented here is realized. In addition, the capture of DM particles in the Sun or 
the Earth is also impossible since the DM will not loose its kinetic energy when scattering 
with electrons in the Sun. Therefore we do not expect high energy neutrino signals from the 
Sun or the Earth either. 
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